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ABSTRACT 

The growth of electrostatic waves near the plasma frequen- 
cy f p due to an unstable electron beam is investigated by 
solving the unmagnetized electrostatic dispersion equation 
numerically. These numerical solutions are compared with 
analytic theories for reactive (or fluid-type) and kinetic ver- 
sions of the beam instability, and for the O’Neil/Malmberg 
connection of the beam and Langmuir modes. Conditions 
for growth significantly above or below f p are given. Three 
general results given are found: (1) The unstable waves do 
not grow on a mode with Langmuir dispersion except in the 
limit of a very dilute beam with growth on O’Neil/Malm- 
berg’s connected mode. (2) The properties of the unstable 
mode depend strongly on beam parameters such as beam 
density, speed and temperature. (3) The frequency of max- 
imum growth frequently lies significantly above or below / p , 
and differs significantly from that predicted by the Lang- 
muir dispersion relation. These results imply important 
consequences for theories of strong turbulence and non- 
linear wave- wave processes, and observational identification 
of f p from observed wave frequencies. 

1. INTRODUCTION 

Electrostatic waves excited near the electron plasma fre- 
quency fp by electron beams axe of widespread interest 
in plasma physics; specific applications include planetary 
foreshocks (Ref. 1), strong Langmuir turbulence (Ref. 2) 
and pulsar magnetospheres (papers this volume). Theoret- 
ical concepts for these electron “beam instabilities” have 
changed relatively little since their initial development 
(Refs. 3-5), in part because most research effort has been 
concentrated on the non-lineax and strong turbulence evo- 
lution of the waves. This paper is, however, specifically 
concerned with the linear properties of the waves and the 
instabilities themselves. 

The two primary characterizations of the beam instability 
are due to Bohm and Gross (Ref. 3) and Briggs (Ref. 4), 
and to O’Neil and Malmberg (Ref. 5). This first character- 
ization corresponds to the physical growth mechanism of 
the waves: viz. the reactive (or non-resonant or fluid-like) 
and kinetic (or resonant or resistive) instabilities in which 
growth occurs due to particle bunching or by inverse Lan- 
dau damping, respectively. These different growth mecha- 
nisms imply different wave spectra for the two instabilities 
(e.g, Ref. 6). The second characterization (Ref. 5) corre- 
sponds to the dispersive character of the growing mode: the 


waves might be “beam modes” with u r ~ it or might lie 
on a mode with approximately Langmuir dispersion formed 
after connection of the ordinary Langmuir mode and the 
beam mode. Although analytic theory (Ref. 7) implies 
that the reactive instability occurs on a beam mode, the 
detailed connection between these two characterizations of 
the beam instability is not clear. One aim of this paper 
is to clarify the connection between these two character- 
izations of the beam instability. This is done by solving 
the exact electrostatic unmagnetized dispersion equation 
numerically, and comparing the numerical solutions with 
analytic solutions. It is shown that the reactive instability 
does tend to occur on the beam mode, but that the insta- 
bility becomes kinetic on a modified extension of the beam 
mode (the ‘modified- beam mode’) long before the O’Neil 
and Malmberg connection between the Langmuir mode and 
the beam mode takes place. The unstable mode is shown to 
be intrinsically beam- associated (i.e., dependent on beam 
parameters, at least) unless the beam is extremely dilute, 
i.e., relative beam densities n*/n o <;10 -4 ; that is, growth 
does not usually occur on a normal mode of the plasma 
in the absence of the beam, such as the Langmuir mode. 
Consequently, there is no a priori reason for growth to oc- 
cur above f p at Langmuir-like frequencies. Conditions for 
reactive or kinetic growth well above and well below / p , as 
observed in the Earth’s foreshock (Refs. 8, 9) are therefore 
investigated here. 

Previous numerical investigations of this dispersion equa- 
tion (e.g., Ref. 10) have paid relatively little attention to 
the dispersive characteristics of the growing waves. Yet the 
detailed dispersion relation of the waves is vital in theo- 
retical descriptions of non-linear wave processes and strong 
turbulence phenomena (e.g., Ref. 2): differences between 
the actual dispersion relation and the assumed dispersion 
relation imply differences in the governing equations for the 
wave fields and possible kinematic difficulties for the wave 
processes, as discussed in Section 4 below. This paper pro- 
ceeds as follows: Section 2 contains a brief description of 
previous analytic theory, while Section 3 contains the re- 
sults of the paper. A more complete version of this work is 
described in Ref. 6. 

2. RELEVANT ANALYTIC THEORY 

Given the longitudinal part of the unmagnetized, spatially 
homogeneous dielectric tensor K L = l+/G+AT e +Aft, where 
the subscripts, i, e, and b refer to the contributions of ions, 
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background electrons and beam electrons, respectively, the 
dispersion equation is K L — 0. Assuming that all species 
have Maxwellian distributions, K a is of the form 


K* 


k 2 V* 


(*( 


U - k-V<, 
V2kV a 



( 1 ) 


where $ is related to the Friede- Conte function, and u/’ po , V a 
and are the angular plasma frequency, thermal speed 
and drift velocity of species a, respectively. Writing u = 
u> r + *7, a wave is said to be resonant with species a if 
= («r-&J’a)/'/2 k V a 1 . The beam instability is said 
to be resonant if < 2, corresponding to the distribution 
function of species a at the phase velocity of the wave being 
2% of its value at v + — v a . 


The reactive instability is derived by ignoring Im K a (i.e., 
Im K a = 0), whence the dispersion equation is a real 
equation for a complex variable u>, thereby implying the 
presence of complex conjugate solutions (Ref. 7). Ignor- 
ing Im K a corresponds to ignoring Landau (and inverse 
Landau) damping; growth is instead due to particle bunch- 
ing. Assuming |£&| 1 (non- resonant) one finds complex 

conjugate roots with maximum growth with u) p = k- 



7r = 


V3 rib 
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The kinetic instability is usually derived by ignoring 
Re(K b , Ki), assuming |£&| < | (resonant), and taking 


7 = 


Im K b 


( 3 ) 


with UJ r determined by the equation Re K t = 0. Growth 
is then due intrinsically to inverse Landau damping, with 
maximum growth rate 


corresponding to the reactive instability are clearly visi- 
ble in Figure 1(a), together with their beam dispersion 
uv ~ kv b and non-resonant character. It is clear that 
the reactive/kinetic transition takes place long before the 
O’Neil and Malmberg (Ref. 5) connection of the Langmuir 
and beam mode takes place in Figure 1(f); furthermore, 
the instability does not take place on the pure Langmuir 
mode, or on a mode with Langmuir-like dispersion until 
greater beam temperatures than those in Figure 1(f) are 
reached. The following changes occur as T* increases in 
Figures 1(b) through (d): (1) the complex conjugacy of 
the beam roots disappears, (2) the dispersion curve of the 
growing mode consists of a portion with beam-type disper- 
sion u) r ~ k V b for wavenumbers smaller than some k m and 
slowly- varying dispersion u; r ~ A for k > A;* , where A is of 
order u p . (3) The wavenumber k m and real frequency u > m 
corresponding to maximum growth increase with Tj,, corre- 
sponding to the peak in the growth rate curve (as a function 
of wavenumber) moving continuously towards and along the 
tj r ~ A portion of the dispersion curve. The instability is 
marginally resonant in Figure 1(b) and strongly resonant 
by Figure 1(c), corresponding to a kinetic instability. By 
Figure 1(d) the portion of the growing mode with beam 
dispersion (i.e., small wavenumbers) is heavily damped, as 
is the Langmuir mode for higher wavenumbers. In contrast, 
the Langmuir mode at small wavenumbers, and the modi- 
fied beam mode at larger wavenumbers are lightly damped 
or growing. As T b increases further, the two lightly damped 
and two heavily damped portions of these modes connect, 
leading to O’Neil and Malmberg’s ‘Langmuir- beam’ mode 
(Gary’s terminology) upon which the instability is found, 
and a heavily damped mode of no further interest. This 
connection is shown in Figures 1(e) and (f). For complete- 
ness note that a progression of figures analogous to Fig- 
ure 1 could also be constructed by either decreasing the 
beam density or decreasing the beam speed, keeping other 
parameters constant. 
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where e = exp(l), on the Langmuir mode. A commonly 
used criterion (e.g., Ref. 7) for determining whether the 
instability is reactive or kinetic is 

** = (—)* 57 >1 ( 5 ) 

n 0 Vb 

or P < 1, respectively. A more accurate empirical criterion 
is given below. 

O’Neil and Malmberg’s (Ref. 5) criterion for whether the 
instability is on the beam-mode or on the connected 
Langmuir- be am mode is 

5 = 2 >/P < 1.47 ( 6 ) 

or a > 1.47, respectively. Gary’s results (Ref. 10) and the 
results below verify this criterion empirically. 

3. SUMMARY OF RESULTS 

3.1 Reactive/kinetic. beam /Langmuir-beam, and disper-_ 

sion relations 


The strong modification of the beam root’s dispersive char- 
acteristics as the instability becomes resonant (and so ki- 
netic rather than reactive) necessitates a clarification of 
terminology: (1) Portions of dispersion relations with o; r ~ 
kv Q are called ‘beam’ type. (2) Due to its origin with the 
beam the strongly modified mode (Figures l(b)-(e)), on 
which growth occurs following modification of the complex 
conjugate beam mode, is called the “modified- beam” mode. 
(3) The connected O’Neil and Malmberg mode is called the 
“Langmuir-beam” mode following Gary (Ref. 10). With 
this terminology the instability is resonant (kinetic) on the 
modified-beam mode and Langmuir-beam mode, and non- 
resonant (reactive) on the beam mode. These identifica- 
tions of the reactive and kinetic instabilities are consistent 
with the variations in the maximum growth rate 7 m as a 
function of relative beam density nb/n 0 in equations (2) 
and (4). Finally, values of P corresponding to the reactive 
instability are P^6, while the kinetic instability is strongly 
resonant when P<> 2. Reactive growth occurs on a beam 
mode which strongly obeys beam-type dispersion. Quali- 
tatively, the modified-beam mode’s dispersion curve con- 
sists of a beam-type portion at low wavenumbers and a 
portion for which u r varies slowly; it is inappropriate to 
approximate the modified-beam mode with the Langmuir 
dispersion relation, particularly as the maximum growth 
rate usually occurs below f p . 


Figure 1 shows the evolution of the dispersion curves for 
a beam with fixed relative density ni/n 0 = 10 -3 and drift 
speed Vb/V e = 10 as the beam temperature Tb/T e is var- 
ied (Tt/Ti = 3 is fixed). The two complex conjugate roots 


The dispersion relation, near where the growth rate is max- 
imum, of the Langmuir-beam mode is also often consider- 
ably different from the Langmuir dispersion relation (Fig- 
ures 1(f) and 2). Figure 2 shows the dispersion relation near 
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maximum growth for beams with various beam densities 
n b /n 0 but constant beam speed v b /V e = 10 and tempera- 
ture T b /T e = 1; the Langmuir-beam connection takes place 
near n b /n 0 = 5 x 10 “ 4 . Stringent conditions on n b /n 0 are 
apparent for the dispersion relation to be closely Langmuir- 
like: n b /n Q should be at least a factor of 10 less than that 
beam density satisfying the O’Neil and Malmberg criterion 
(6), i.e., n b /n 0 < 5 x 10" 5 . This rule should be of great 
practical importance: for the above warm, relatively small 
velocity beam one requires n b /n Q < 5 x 10“ 5 , and for faster, 
colder beams one requires even more dilute beams; however, 
in most laboratory and space science applications (e.g., the 
Earth’s foreshock, but not Type III solar radio bursts) one 
expects n b /n Q ^ 10“ 3 . 

This sub-section may be summarized as follows: (1) The 
unstable waves do not have the Langmuir dispersion rela- 
tion, except in the limit of a very dilute beam with growth 
on the Langmuir-beam mode: n t /n 0 <5x 10~ 5 for v b /V e 
^10 and T b /T e ^ 1. (2) In most applications growth should 
occur on the beam mode (reactive) or modified- beam mode 
(kinetic). (3) The properties of the unstable waves depend 
strongly on the beam parameters, such as beam density, 
speed, and temperature. (4) The frequency of maximum 
growth frequently lies either significantly below or above f p , 
and is significantly different from the frequency predicted 
by the Langmuir dispersion relation. 

3.2 Growth significantly above and below f r 

Figures 1 and 2 show that the frequency of maximum 
growth for the beam instability often lies below / p , despite 
the conventional belief (due to the idea that growth oc- 
curs on the unmodified Langmuir mode) that growth oc- 
curs above / p . Effective growth well above or below f p 
is due to the strong dependence on beam parameters of 
the dispersion relations of the beam, modified-beam and 
Langmuir-beam (unless the beam is very dilute) modes, 
and the non- Langmuir nature of these modes’ dispersion 
relations. Figure 3 illustrates the roles of beam speed, tem- 
perature and density in determining the real frequency at 
maximum growth a ; m (Ref. 9). Growth significantly be- 
low oj p generally (except as in Figure 2) requires growth 
to occur on the beam (reactive) or modified-beam (kinetic) 
modes. For high beam speeds ( v b /V t > 20) u m becomes 
approximately constant, equal to u> n say, and dependent 
only on n b /n a . For fixed beam density and temperature, 
decreasing the beam speed causes u m to increase above u> n 
once Ufr/V’e-^IO, and subsequently to decrease towards zero 
frequency once v b /V e <&2, provided the beam temperature is 
sufficiently small. This increase in above co„ increases 
with decreasing n b /n 0 until n b /n 0 ~ 10~ 5 . Increasing the 
beam density corresponds to moving the entire curve in 
Figure 3, especially u; n , to lower frequencies: note that 

for n b /n 0 ^ 10~ 3 ,u? n «£0.99u; p . For warm beams (e.g., the 
dashed curves) damping overcomes growth before u: m can 
increase substantially above u n , and certainly before u> m 
decreases towards zero frequency. This Figure provides a 
natural explanation for waves observed at frequencies sig- 
nificantly above and below u ) p in the Earth’s foreshock (see 
Refs. 8, 9). 

There are therefore two primary ways to obtain growth sig- 
nificantly below u? p : dense beams or cold, slow, relatively 
dilute beams. Growth significantly above u p requires rela- 
tively slow (2^Ufc/V e ^10) and dilute (since otherwise u> n is 
well below v p ) beams which may be either cold or warm. 
Detailed conditions for reactive or kinetic growth well above 
or below u p are given in Ref. 9. 


4. DISCUSSION 

Due to length restrictions the discussion is restricted to 
the following two points: (A) The results of this pa- 

per indicate that, due to the assumption of the Langmuir 
dispersion relation in derivations of the Zakharov equa- 
tions (e.g., Ref. 2), and the neglect of particle bunching 
effects therein, the Zakharov equations are valid only un- 
der stringent circumstances: n b /n 0 <£ 10“ 4 for v b /V e > 10 
and T b /T e % 10. (B) The generally non-Langmuir disper- 
sion relations, and tendancy for maximum growth to oc- 
cur below u> p , of waves growing by the beam instability 
imply kinematic difficulties for non-linear wave processes. 
For example, frequency and wave vector conservation for the 
Langmuir wave decay L -> V + S (L denotes the electron 
.plasma wave, while S is an ion acoustic wave) is impossi- 
ble when nft/n o ^10 -3 since u) m is below u ; p , yet a; Jr, (the 
backward propagating wave) is greater than u p (Ref. 6). A 
more detailed discussion is given in Ref. 6. 
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Figure 1. Plots of the dispersion relation and growth rate for the growing mode (full line), 
damped beam mode (dashed line), Langmuir mode (dotted line) and Langmuir-beam 
mode (dash-dot line) for various beam temperatures. For parts (a) to (f ) T b /T e equals 
0.017, 0.067, 0.33, 1.0, 1.5, and 1.67, respectively. The fixed beam parameters are 
n b ln 0 = 10 -3 , v b /V t = 10, and T e /Ti = 3. 
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Figure 2. The dispersion relation of the Langmuir-beam 
mode as a function of n b /n 0 : 5 X 10 -4 (full 

line), 10 -4 (dashed), 5 x 10 -5 (dash-dot), and 
10 -5 (dotted). Here v b /V e — 10, T b /T e — 1, 
and T e /T, = 3. Connection of the Langmuir 
and beam modes occurs near n b /n 0 = 5 x 10“ 4 . 
The mode has the dispersive characteristics of 
the Langmuir mode for rib/n 0 < 10 -5 . 
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Figure 3. The frequency of maximum growth, w m , as a 
function of v b /V e for n b /n 0 = 10 -2 and various 
beam temperatures. For high beam speeds u m 
is approximately constant, equal to o; n , and pri- 
marily dependent only on the beam density. If 
the beam is sufficiently cold, u> m first increases 
above and then decreases below u n as v b /V e de- 
creases. 
















